How natural is a small but nonzero cosmological constant? 
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Abstract 

Based on our previous attempt, we propose a better way to understand 
a small but nonzero cosmological constant, as indicated by a number of 
recent observational studies. We re-examine the assumptions of our model 
of two scalar fields, trying to explain the basic mechanism resulting in a 
series of mini-infiations occuring nearly periodically with respect to Int 
with t the cosmic time. We also discuss how likely the solution of this type 
would be, depending on the choice of the parameters. 

A growing number of different observations, notably the recent determination of 
the Hubble constant Hq [0, seem to point to a suggestion that there is a pos- 
sible small but nonzero cosmological constant A 0, with Q/^ = A/pcr~l; where 
Per = (3/87rG)ifQ. This may not, however, be readily acceptable from a theo- 
retical point of view, because introducing A has been considered to be highly ad 
hoc. Contrary to this long-held prejudice, on the other hand, it is widely recog- 
nized that a cosmological constant is an indispensable ingredient in many of the 
theoretical models of unification.]] Unfortunately, they tend to predict A larger 
than the observed value, or its upper bound, by as much as 120 or so orders 
of magnitude. One of the possible ways out is to devise a theory in which the 
cosmological constant is not a true constant but decays like ~ t~^, with t the 
cosmic time @,|]. 

Notice that the theoretically natural size of A is of the order one in the Planck- 
ian unit system with SnG = 1,| while the present age of the Universe to ~ 10^°y 
being of the order of 10^°. In this scenario of "a decaying cosmological constant," 
today's cosmological "constant" is small <jlO~'^^^ only because our Universe is 
old. No unnaturally extreme fine-tuning of parameters is called for. 

The scenario has been often formulated based on the models in which a scalar 
field plays a role; the time-dependent effective cosmological constant Aefj(t) is in 
fact the energy density of this scalar field that couples to the ordinary matter 
only as weakly as gravity. 

*E-niail address: ysfujii@tansei.cc.u-tokyo.ac.jp 

^Some authors ||^ assume A = in the starting Lagrangian. This is protected, however, by 
no known symmetries against any likely perturbations. 

■'■By also choosing 7i = c = 1, the units of length, time and the energy are 8.09xl0^'^^cm, 
2.70xl0~"*3sec, and 2.44x lO^^GeV, respectively. 



This model implies, however, a complete absence of the cosmological constant. 
What is indicated by the observation is, on the contrary, the presence of a flat 
portion in the energy density as a function of t, a deviation from a smooth falling- 
off ~ We came across, however, to a model which, by introducing another 
scalar field, would result in occasional fiattenings of Aeff(t) thanks to a nonlinear 
nature of the cosmological equations ||^ . The purpose of this note is to provide a 
simpler understanding of the mechanism. 

It seems appropriate here to state our attitude. In view of the lack of the 
complete theoretical framework to derive all the details of the final results, we 
follow a heuristic approach, trying to see what the effective theory in 4 dimensions 
should be like, if it is to fit to what appears to be the effect of a small but nonzero 
cosmological constant. As it turns out, this is highly nontrivial, if the model 
is somehow related to modern unification theories. We list some of the main 
assumptions in Ref. P]. 

First we assume the presence of a scalar field of the dilaton-tjpe, having a 
non-minimal coupling, which is chosen, for the sake of simplicity, to be a Brans- 
Dicke type; (fP'R.^ We then apply a conformal transformation (Weyl rescaling) 
to remove the non-minimal coupling. We do this for the technical convenience, 
at the moment, though the correct conformal frame (CF) should be selected 
according to what clock we use to describe the evolution of the Universe. In 
this connection we should notice that none of the realistic theories of gravity is 
conformally invariant, and that our conclusion remains true also in the original 
CF in which the non- minimal coupling is present. 

As an important consequence of this transformation the A term in the original 
CF is converted to a potential of (j) of the type Ae~°"/'' where k is a constant^] 
while 0" is a transformed scalar field appropriate in the new CF. The a field rolls 
down the slope toward infinity, ensuring the Aeg, essentially the energy of a, to 
fall off like ~ after the inflationary era. 

As a next step we introduce another scalar fleld $ which has a speciflc in- 
teraction with cr but couples to conventional matter flelds also as weakly as the 
gravitational interaction. We discovered an example of the interaction such that 
Aeff, which is now the total energy density of the 0"-$ system, shows a repeated 
occurrence of leveling-off superimposed on the overall smooth fall-off ~ com- 
ing barely to exceed the normal matter density, hence acting as a cosmological 
constant. In accordance with this the scale factor a{t) exhibits an extra accel- 
eration deviating from the overall smooth behavior ~ t^/^ or ~ t^f^. A typical 
solution of the cosmological equations was shown in Fig. 3 of Ref. . 

The basic cosmological equations are (/c = 0) 

= + Pm, (1) 



^Our is, following the standard notation in the conventional relativistic field theory, related 
to the original notation Lp in Ref. 0| hy ip — cjp' / diuj . 

^See Ref. on how k is constrained in order for the results to be consistent with realistic 
cosmology. 
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a + 3H& — exp(— ct/k) 



A 1 fU dU\ 

— + -m^$^ — 

K 2 \k da I 



^ + 3iJ<i> + exp(-or//t)m^^7$ = 0, 



where 



and 



with 



2 2 



V{(^, a) = exp(-cr/K) 



(2) 



(3) 



(4) 



(5) 



(6) 



U{a) = 1 + Bsm{uja). 

Here K,m,B and u are constants. The exponential factor exp(— ct/k) (~ 
comes typically from the Weyl rescaling, thus transforming the A-term into the 
potential, as mentioned before. The sinusoidal dependence in @, on the other 
hand, has been introduced on the try-and-error basis. We later discuss how this 
specific form is favored. It is crucial to assume that the conformal transformation 
property of $ is such that the same factor exp(— ct/k) appears in front of the 
$^f/(cr) term. 

For the matter density pm we assume the mixture of relativistic and non- 
relativistic matters: 

Pm = Prfl"^ + Pnra"^- (7) 

We admit that the result depends heavily on these assumptions. In view of 
the huge discrepancy of 120 orders of magnitude between conventional theory 
and the observation, however, variety of models seem to deserve consideration as 
working hypotheses. Notice that all of our parameters are essentially of the order 
one in Planckian units, appealing to theoretical naturalness. 

One might still argue that we are introducing as much as what we want to 
come up with. We emphasize, however, that it is far from trivial to make a right 
choice on what to be introduced; otherwise it will not make sense no matter how 
much we bring in. 

On using the new time variable r defined by 



and also defining 6(r) by 

eqs. (|1])-® ai'e put into 

-a' ^ 4 
2 2 



r = Int, 



a = e 



3b'^ = la'^ + U'^ + t''{V + p, 



a" + (36' -l)a' -t^ exp(-(T//t) 



A 1 2^2 (U dU\ 

— + -m^<l>^ — 

K 2 \K da J 



0, (9) 



$" + (36' - 1)$' + t^ exp(-a/K)m2f/$ = 0, 



(10) 
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where ' means differentiation with respect to r. Notice the exphcit occurrence of 
the time variable on the right-hand sides. 

It is also worth noticing that (^-(p!OD allow the asymptotic solution for t ^ oo 
if U is chosen to be constant, namely B = 0: 



a{t) = 



where A is an integration constant while 



or 6(r) = -r, (11) 

(12) 
(13) 




<|)(r) = Ae ■^/^sin(mr), 



12^?^? 1 

We learn that the $ field, if decoupled from cr, would play no role in the 
asymptotic era and that r might be a useful variable. From the solution (|1^) also 
follows that the combination 

F{t,o) = t^exp(-a/fi;) 

tends to a constant (= 2/t^/A) if [/ = 1. This implies that with non-constant 
f/(cr) this combination might be nontrivial. We show that this is indeed the case. 

A typical solution obtained numerically is shown in Fig. 1, another example 
with parameters somewhat different from those used in Fig. 3 in Ref. [^]. Notice 
that we chose the "initial time" ti = 10^°, because, though the real classical 
cosmology had begun much earlier, we can conveniently avoid more details on 
the inflation era and the ensuing reheating process. 

We first find in the plot (a), as already alluded, the scale factor a{t) shows a 
series of "mini-inflations," each implying a rapid increase lasting during an opoch 
which is "short" in terms of r, but could be quite "long" if it is measured in 
the ordinary time t; nearly as comparable as t itself. One of such epochs has 
been chosen to include the present time with A = logt ~60, corresponding to 
t ^ 10^°y. (The scale factor resumes a usual expansion immediately beyond the 
frame.) 

In the plot (b) we notice that ps, the total energy of the 0"-$ system, and 
Pm; the ordinary matter energy density, fall off like ~ as an over- all behavior, 
but interwinding each other. A closer look reveals that a mini-inflation occurs 
whenever exceeds pm- 

Fig. 2 is the same plot as Fig. 1 (b) presented in a magnified scale around 
the present time; showing that ps surpasses pm, remaining nearly constant for a 
while, hence imitating a cosmological constant with the size basically of the same 
order of magnitude as ~ pcr- 



We also notice in Fig. 1(a) that each of these anomalous behaviors takes place 
toward the end of the "dormant period," during which both of the scalar fields 
come almost to standstill. To understand this behavior we first point out that 
the dormant period and its repetition are primarily due to the dynamics of the 
(7-$ system; the "back-reaction" from the cosmological expansion has a rather 
minor effect. 

Fig. 3 shows an example in which the cosmological part is cut off with the 
same parameters but with 36' — 1 replaced by 0.5. In spite of some differences, 
which, representing how much the cosmological effect could be, will be discussed 
shortly, comparing Fig. 1(a) and Fig. 3 is sufficient to convince ourselves that 
the "recycled dormant periods" could take place even without cosmological effect. 
From this point of view, we now focus upon more detailed analysis of the solution 
in the isolated a -$ system as a simplified mathematical model. 

The initial value cxi = 6.75442 implies that uai = 27r x 10.75; the a field 
starts at one of the potential minima as given by sin(a;cr) if viewed in the a 
direction (see Fig. 2 of Ref. [§), but on the slope in the direction of $. Also the 
"initial time" chosen to be ti = 10^° implies 2ti = 46.052, and hence F{ti, ai) = 
exp(— (Ti/k) = exp(2ri — cti/k) = exp(3.303) = 27.19, which is quite large. As 
a result, a is pushed forward strongly. In this sense the system started from a 
"catapulting stage." The rapid increase of a, however, makes F{t, a) small, as 
will be found by comparing the curve of a and the straight line 2kt, also shown in 
Fig. 3. Soon a is nearly free, going further until it is decelerated by the frictional 
term 36' — 1 ~ 0.5, finally to be trapped to another minimum of sm{uja). 

On the other hand, the $ field, having pearched on the middle of the potential 
slope, is also catapulted downward, passing the central valley $ = past, until the 
force dwindling again due to the decrease of F[t, a) and the cosmological friction 
stop it to an almost complete halt, hence the dormant period. The energy density 
Ps still continues to decrease according to ~ before it stays constant.]^ 

Now with a virtually unchanging a, the increasing r makes F{t, a) non- 
negligible again, bringing the system back to the catapulting stage from which we 
started before. In this way the dormant period may repeat itself nearly periodi- 
cally with respect to r (instead of t) , if the field configurations match suffciently 
close to the previous values. 

The "recycling," however, may fail if the matching turns out incomplete. Fig. 
3 is in fact one of the patterns of such "short" recycling encountered most com- 
monly. Toward the end of a dormant period, a is "released" off the track before 
it is kicked hard by the force which has been building up. The system enters 
into an asymptotic behavior in which a grows linearly while $ decreases slowly 
toward $ = 0. If the same behavior occurs in the cosmological setting, the scale 
factor increases smoothly, resulting in no effective cosmological constant today. 

On the other hand, much "longer" recycling is also rather common, as shown 
in Fig. 4. Suppose a recycling process lasts suffciently long in the cosmological 
system. Then, as in the isolated 0"-$ system, would stay nearly constant 

II See also Ref. |^ for a similar behavior. 
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toward the end of the dormant period, surpassing pm, hence playing the role of a 
cosmological constant. When the scalar fields start moving as the factor F[t, a) 
increases, however, ps begins to decrease, eventually nosediving beneath p^. The 
Universe resumes an ordinary expansion again. This explains the behaviors shown 
in Figs. 1-2. 

A question then arises how likely the solutions of sufficiently long recycling 
could be. An idea on the answer may be obtained again by studying the isolated 
cr-<l> model. We surveyed solutions of the isolated 0"-$ system by changing one 
of the initial values, $i for example, keeping other parameters and initial values 
fixed.Q Solutions with shorter recycling exhibit basically the same patterns as in 
Fig. 3, while Fig. 4 is an example of sufficiently long recycling. Combining the 
solutions, we plotted in Fig. 5, the time of the end of recycling, te, against $i 
varied discretely. In spite of apparently rampant variation, we obtain solutions 
of long recycling (Ae > 62) for 8 out of 23 choices of $1. This, together with 
other limited but similar examples, seems to be an encouraging sign that the 
occurrence of the continued recycling is reasonably likely. The same "optimistic" 
view applies also if the cosmological effects are fully included. 

We emphasize that the presence of minima of the potential as a funciton of a 
is crucial. The leveling-off behavior is triggered by trapping a. The form sin(co'cr) 
is favored because it is ready to trap cr virtually at any time. Any other potential 
will probably be acceptable, from a "phenomenological" point of view, if it shares 
this property. 

We have shown that a small but nonzero cosmological constant as required 
by the observations could result rather naturally due to a nonlinear nature of 
the scalar field equations. As a favored coupling, we tentativley suggested a 
Sine-Gordon-like interaction, though its real origin is yet to be discussed 0. 
Furthermore having introduced two scalar fields, we have too many parameters, 
including the initial conditions, to allow unique predictions, or even a systematic 
survey of the solutions. In this sense our conclusion is still preliminary. The 
example in Figs. 1-2 merely illustrates how the results can be realistic, leaving 
many details yet to be worked out. We nevertheless have promising indications 
that the desired result comes about quite likely. 

We certainly have to adjust some of the parameters in order to bring a mini- 
infaltion to the epoch including the present time, for example. The extent to 
which we are supposed to fine-tune them is rather mild, however In other words, 
detailed analyses of various cosmological parameters at present, like HQ,tQ,Q\ 
and go 5 will serve to constrain the parameters of the theory, as will be discussed 
elsewhere. 

As one of the generic consequences of the present mechanism, we anticipate 
some past epochs to have emerged with significant amount of In the solution 

**Presentation of the result has been made simpler by choosing parameters different from 
those used in Figs. 1-3, but the same as in Fig. 4 except for $1. 

ttChanging $1 from 0.212 to 0.210 in the solution of Figs. 1-2 would shift f^A = 0.67 at 
A = 60.15 {t = 1.21 X W^^y) to 0.71. However, $1 = 0.200 gives I^a - 1 - lO^^o q^^ ^j^g ^^y^^^ 
hand, $1 = 0.2115 yields Ha = 0.73, Hq =79 km/sec/Mpc at t = 1.5 x lO^^y 
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of Figs. 1-2, for example, we find JIa » 1 for 27<,A <,39, while we have avoided the 
same at the epoch of primordial nucleo-synthesis (A ~ 45). This may illustrate 
how the parameters and the initial conditions can be determined, in principle, also 
by looking into details of the past cosmological histories, which should deserve 
future studies. 

We confined ourselves to the "primordial" cosmological constant prepared 
in the starting Lagrangian. It is yet to be seen if the same mechanism applies 
successfully to the vacuum energies associated with cosmological phase transitions 
at later times. 
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Figure 1: An example of the solution of (8)-(10). (a) Upper plot: 6 = Ina 
(solid), a (dotted) and 2$ (broken) are plotted against A = logi = 0.434t. 
The present age of the Universe supposed to be (1.0- 1.5) xlO^°y corresponds to 
60.0-60.2 of A in units of the Planck time. The parameters were chosen to be 
A = 1, «; = 0.158, m = 4.75,5 = 0.8, a; = 10 in the Planckian units. The initial 
values chosen conveniently at t\ — 10^° are a — = 6.75442, (j\ = 0, $i = 
0.212, $1 = 0, pri = 2.04 X 10-2\ pnri = 4.46 X 10-^^; the last two being adjusted 
to give the "equal time" Aeq ~ 55. The value of a\ corresponds to starting at a 
minimum of sin(cL;a). (b) Lower plot: ps (solid), the total energy density of a and 
$, and pm (dotted), the matter energy density, against A = logt. 
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Figure 2: The same plot as in Fig. 1(b) but in a magnified scale of A around 
the present time. We find Qa — 0.67 and Hq — 81km/sec/Mpc at A = 60.15 
{t = 1.21 X lO^V)- 
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Figure 3: An example of the solution in the isolated (7-$ system, in which re- 
cycling of the dormant periods ends prematually at A ~ 41. a (solid) and 2$ 
(dotted) are shown against A = logt. Also shown is 2kt (broken) to be compared 
with a. All the parameters remain the same as in Figs. 1-2, except for 36' — 1 
replaced by 0.5. 
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Figure 4: An example of the solution in the isolated a-^ system, showing long 
recycling. We choose m = 5.0 with other parameters as well as the symbols the 
same as in Fig. 3. The initial values at Ai = 10 are cri = 8.0, di = 0, $i = 
2.1, $1 = 0.19, somewhat different from those in Fig. 3. 
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Figure 5: The time Ae for the end of recychng in the isolated (7-$ system is 
plotted against one of the initial values $i, varied with spacing 0.05. The other 
initial values and the parameters are the same as in Fig. 4. The arrows indicate 



Ae > 62. 
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